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\ ABSTRACT. In this paper we apply the semidefinite programming approach de- 

veloped in (2) to obtain new upper bounds for codes in spherical caps. We com- 
pute new upper bounds for the one-sided kissing number in several dimensions 
where we in particular get a new tight bound in dimension 8. Furthermore we 
show how to use the SDP framework to get analytic bounds. 
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Dedicated to Eiichi Bannai in occasion of his 60th birthday 

oj . 1. Introduction 

Let S n ~ 1 denote the unit sphere of the Euclidean space M n . The spherical cap 
with center e G S n and angular radius <fi is the set 

^ ■ Cap(e, (f) = {x € : e • x > cos <fi}. 

^ \ We consider the problem of finding upper bounds of the size of a code C contained 

in Cap(e, 4>) with minimal angular distance 6. Following notations of |3], the max- 
imal size of such a code is denoted by A(n, 6, <f>). Many reasons to consider this 
problem are exposed in |3], e.g. upper bounds for spherical codes can be derived 

q . from upper bounds for spherical cap codes through the following inequality: 

23 : A(n,9) < A(n,e,4>) 



X 



voK^"- 1 ) ~ vol(Cap(e, <j>)) 

where A(n, 8) stands as usual for the maximal size of a spherical code with mini- 
J> \ mal angular distance 9. 

Moreover, it is a challenging problem, because the so-called linear program- 
ming method does not apply to this situation. In coding theory many of the best 
upper bounds are consequences of the so-called linear programming method due 
to P. Delsarte. This method gives upper bounds for codes from the solution of a 
certain linear program. It can be applied to symmetric spaces and has been success- 
fully used to deal with two-point homogeneous spaces like the unit sphere S n ~ l 
(1H, 10, ifm and the survey [7, Chapter 9]), or with symmetric spaces which are 
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not two-point homogeneous like Grassmannian spaces ([1]). However the method 
is not applicable to spaces which are not symmetric spaces like spherical caps. 

In this paper, we show that the approach developed in [2] based on semidefinite 
programming can be applied to the above problem. It turns out that it gives good 
numerical results. In particular we obtain improvements in the determination of 
the so-called one-sided kissing number, corresponding to <fi = tt/2 and 9 = 7r/3, 
and denoted by B(n) after fl4l . 

Let us describe briefly the idea underlying our approach. The isometry group 
of Cap(e, 0) is the group H := Stab(0(M n ), e) stabilizing the point e in 0(M n ). 
This group acts on the space Pol<d(S' n_1 ) of polynomial functions on the unit 
sphere of degree at most d. In the decomposition of this space into irreducible sub- 
spaces some irreducible subspaces occur with multiplicities. To each irreducible 
subspace with multiplicity m we can associate an m x m matrix Y whose coeffi- 
cients are real polynomials in three variables (u, v, t) and have an explicit expres- 
sion in terms of Gegenbauer polynomials. Each matrix Y satisfies the positivity 
property: 

For all finite C C S" n_1 , ^ Y{e ■ c, e ■ c', c • c') h 0, 

(c,e')GC* 2 

where u y 0" stands for "is positive semidefinite". 

We want to point out that one can consider other metric spaces X with isometry 
group in this framework. Only the expression of the matrices Y will depend on 
the specific situation. For a symmetric space X the multiplicities in the irreducible 
decomposition are equal to 1. Hence the matrices Y have size 1 x 1. So we recover 
the classical positivity property of zonal polynomials which underlies the linear 
programming method. 

The paper is organized as follows: Section |2] recalls the needed notations and 
results of [2]. Section [3] states the semidefinite program (SDP for short) which 
obtains an upper bound for A(n, 9, 4>) and presents numerical results. Section 0] 
translates the dual SDP into a statement on three variable polynomials, and contains 
more material on orthogonality relations, positivity property and other classical 
material which might be of independent interest. 

2. Review on the semidefinite zonal matrices 

We start with some notations. The standard inner product of the Euclidean space 
IR n is denoted by x ■ y. The orthogonal group 0(M n ) acts homogeneously on the 
unit sphere 

S 11 - 1 ■= {x G R n : x ■ x = 1}. 

The space of real polynomial functions on 5 n_1 of degree at most d is denoted by 
Pol<d(5' n_1 ). It is endowed with the induced action of 0(M n ), and equipped with 
the standard 0(E n ) -invariant inner product 

(1) (/,<?) = — / f(x)g(x)du; n (x), 

w n JS"- 1 
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where uj n is the surface area of S n 1 for the standard measure duo n . 
It is a classical result that under the action of 0(M ra ) 

(2) Pol^S"" 1 ) = H% JL H? JL . . . ± H2, 

where iJP is isomorphic to the 0(M n ) -irreducible space of homogeneous, har- 
monic polynomials of degree k in n variables, denoted by Harm^. For the dimen- 
sion of these spaces we write hu '■= dim(Harm^). 

For the restricted action of the subgroup H := Stab(e, 0(M n )), introduced 
above, we have the following decomposition into isotypic components: 

(3) Poles' 1 - 1 ) = J J. li J_ . . . JL I d , 
where 

l k ~ (d- fc + lJHaimjp 1 , k = 0,...,d. 
More precisely, l k decomposes as 

(4) I k = Htf±...±H%\ 

where, for i > k, H^J is the unique subspace of Hf isomorphic to Harm^ . 
The following construction associates to each T k a matrix- valued function 

(5) Zl : S"- 1 x 5 n_1 -» 

which is uniquely defined up to congruence. Let (e* l5 2 , • • • ; e fc , n -i) be an 
orthonormal basis of H^ +s . We assume that the basis (e* J 1<i<ft »-i is the image 
of( 

e o,i)i<i<h™- 1 some ^-isomorphism : H kk — > H k k ^ s . Then, define 



/ e o,i(^) ••• e 5 )h n-i(a?) ^ 



^r 1 



and 

(6) ZZ(x,y) := EflWEMy)* e 



-fc+l) x (d-fc+l) 



One can prove that, for all g £ H, Z k (g(x), g(y)) = Z%(x,y). As a conse- 
quence, the coefficients of Z k can be expressed as polynomials in the three vari- 
ables u = e ■ x, v = e ■ y, t = x ■ y. More precisely, let Y k n (u,v,t) be the 
(d — k + 1) x (d — A; + 1) matrix such that 



(7) Z%(x,y) = Y k n (e-x,e-y,x-y). 

1 U., L 

fc 



We denote the zonal polynomials of the unit sphere S n 1 by PP. In other words, 



if n > 3, then P]}(t) is the Gegenbauer polynomial of degree A; with parameter 
n/2 — 1, normalized by the condition P k {l) = 1. If n = 2, then -P^(i) is the 
Chebyshev polynomial of the first kind with degree k. We give in [2, Theorem 3.2] 
the following explicit expressions for the coefficients of the matrices Y k n : 
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Theorem 2.1. We have, for all < i,j < d — k, 



(8) (Y£). .(u,v,t) = \ hj P? +2k { U )P^ 2k {v)Ql-\u,v,tl 



where 

Ql~\u, v, t) : = ((1 - u 2 )(l - v 2 )) k/2 P?- 1 ( - f ~ UV =) , 

and 



^ _ U n (J n+2 fc-l /^ n +2fc^n+2fcNl/2 



We recall the matrix-type positivity property of the matrices Y fc n which underlies 
the semidennite programming method: 

Theorem 2.2. For any finite code C C 5 n_1 , 

(9) Yk(e-c,e-c',c-c')hO. 

(c,c')ec 2 

Proof. We recall the straightforward argument: 

£ ^(c,c') = (E^( c ))(E^( c ))^ - 

(c,c')GC* 2 ceC cGC 

□ 

3. SEMIDEFINITE PROGRAMMING BOUND FOR CODES IN SPHERICAL CAPS 

Let C C Cap(e, 0) be a code of minimal angular distance 0. Define the domains 
A and Ao by 

A := {(u, v, t) : cos 4> < u < v < 1, 

-l < t < cose, 

1 + 2uvt-u 2 - v 2 - t 2 > 0}, 

and 

Ao := {(n, u, 1) : cos 4> < u < 1}. 
The two-point distance distribution of C is the map y : A U Ao — > K given by 

y(u, v, t) = — — cardl (c, c ) € C : e • c = u, e ■ c = v,c ■ c = 1 1, 

card(C7) 

where 

II it u = v. 

We introduce the symmetric matrices Y^(u, v, t) defined by 

Y n k (u, v, t) -.= i v, t) + y fc >, u, t)) . 
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Then, © is equivalent to the semidefinite condition 

^2 y(u,v,t)Yl{u,v,t) ^0. 

(u,v,t)eAuA 

For any d > 0, the y(u, v, i)'s satisfy the following obvious properties: 
y(u, v,t) > for all (u, v,t) G A U A , 
y(u, v, t) = for all but finitely many (u, v, t) € A U Ao, 

2/( u >M) = l, 

(u,w,i)eA 

^2 y(. u ,v,t) = card(C), 
(«,«,t)eAuA 

y(u, u, t>, i) >: for k = 0, . . . d. 

(ti,(j,t)eAuA 

Hence a solution to the following semidefinite program is an upper bound for 

A(n,e,<f>). 

max { 1+ ^2 y( u , v ,t) '■ 
(«,u,t)eA 

y(u, v,t)>0 for all (u, v, t) G A U A , 

y(u, v, t) = for all but finitely many (u, u, t) G A U Ao, 

^2 y(u,u,l) = l, 
(«,M,i)eA 

v, t)Y^(u, v,t) ^ for all k = 0, . . . , d}. 

(«,u,t)eAuA 

As usual, the dual problem is easier to handle. The duality theorem says that any 
feasible solution of the dual problem provides an upper bound for A(n, 9, (/)). For 
expressing the dual problem we use the standard notation (A, B) = Trace(AB*). 

Theorem 3.1. Any feasible solution to the following semidefinite problem provides 
an upper bound on A(n, 9, (p). 

min { 1 + M : 

F k h0 far all k = 0,...,d, 

d 

(10) J2(F k ,Yl{u,u,l)} <M for all (u, u, 1) G Ao, 

fc=0 

d 

Y,{Fk,Yl(u,v,t)) < -1 for all (u,v,t) G A} 

fc=0 



In order to make use of this theorem in computations we follow the same line as 
in (21 Section 5]. A theorem of M. Putinar ( ifTTl ) shows that the two last conditions 
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Table 1. Bounds on B{n). 



can be replaced by: 

d 

^2( F k,Yl(u,u,l) = M - q (u) -p(u)qi(u) 

d 4 

^2( F k,Yk(u, v,t)) = -1 - r (u,v,t) - ^2pi(u,v,t)ri(u,v,t) 

k=0 1=1 

wherep(n) = — (u— cos <p){u— 1), p\ =p(u),p2 = p{v),ps = — (t+l)(t— cos 9), 
Pi = —{v? + v 2 + t 2 ) + 2uvt + 1, and the polynomials Qi(u), < i < 1 and 
ri(u,v,t), < i < 4 are sums of squares of polynomials. If we set the degree 
of those polynomials to be less than a given value N, and fix the parameter d, we 
relax (fTOl to a finite semidefinite program. 

In the most interesting case cos 4> = and cos 8 = 1/2, corresponding to the so- 
called one-sided kissing number B(n), we obtain the computational results given 
in Table 1. For our computations we chose the parameter d = N = 10. 

In this table, the values in the column of the best lower bounds known cor- 
respond to the number of points in an hemisphere from the best known kissing 
configurations, given by the root systems D%, D4, D§, E§, E-j, Eg. 

Our method gives a tight upper bound in three cases. In dimension 3 we get 
with parameters d = N = 4 the bound B(3) < 9.6685 and hence we recover the 
exact values B(3) = 9 first proved by G. Fejes Toth ( ifTUll ). In dimension 4 we get 
with parameters d = N = 6 the bound B(4) < 18.5085 and hence we recover the 
exact value B(A) = 18 first proved by O.R. Musin (O). In dimension 8 we find a 
new tight upper bound. The famous configuration of 240 points of S 7 given by the 
root system Eg is well known to be an optimal spherical code of minimal angular 
distance 7r/3, which is moreover unique up to isometry. Optimality is due to A.M. 
Odlyzko and N.J.A. Sloane ( lfl6lD . and independently to V.I. Levenshtein ([13]), 
uniqueness is due to E. Bannai and N.J.A. Sloane ([6]). From these 240 points we 
get a code of the hemisphere as follows: Take e among these points, then the subset 
of those points lying in the hemisphere with center e consists in 183 points. We 
obtain a bound of 183.012 with d = N = 8 in our computation. Hence, it proves 
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that it is a maximal code of the hemisphere, in other words that 

B(8) = 183. 

It is reasonable to believe that the configuration of 183 points of Eg is unique up to 
isometry. Unfortunately we were not able to prove it. 

4. Polynomials 

4.1. Polynomial restatement of the SDP bound for codes in spherical caps. 

We want to give an equivalent expression of the bound provided by Theorem B.ll in 
terms of polynomials. Such an expression will be useful to prove analytic bounds 
without the use of software for solving semidefinite programs, just like in the case 
of the linear programming (LP) bound (see e.g. |[T6l ). Moreover, we aim at setting 
bounds in the form of explicit functions of cos 9 and cos <f). We start with a lemma 
which shows that any polynomial in the variables u, v, t can be expressed in terms 
of the matrix coefficients of the Y^(u, v, t). In our situation it suffices to restrict to 
polynomials which are symmetric in u, v. We introduce the following notation: 

R d := {F e R[u,v,t] : F(u,v,t) = F(v,u,t),deg (u>t) (F) < d,deg t (F) < d}, 

where deg( u t ) stands for the total degree in the variables u, t. 

Lemma 4.1. Let F(u,v,t) € Rd- There exists a unique sequence of d + 1 real 
symmetric matrices (Fq, F±, . . . , Fj) such that F^ is a (d — k + 1) x (d — k + 1) 
matrix and 

d 

(11) F(u,v,t) = J2( F k,Yl(u,v,t)}. 

k=0 

We shall say that (Fq, . . . , F^j are the matrix coefficients of F. 

Proof. The polynomials v, t) have degree k in the variable t. Hence, 

F(u, v, t) has a unique expression of the form 

d 

F(u,v,t) = ^2q k (u,v)Ql~ 1 (u,v,t), 

k=0 

where qk(u,v) is symmetric in u,v and has degree in u at most d — k. Since 
P? +2k {u) has de gree i, q k has a unique expression as a linear combination of the 
products \ i jPl l+2k (u)P™ +2k (v) for < i, j <d - k. Thus, there is a symmetric 
{d — k + 1) X (d — k + 1) matrix F k so that 

q k (u,v)= Y, (Fkk 3 K J Pr +2k ^)P] l+2k (v). 
0<i,j<d-k 

Since one can write Y k (u,v,t) as Q™- 1 (u,v,t)(\ i jP* l+2k (u)P! j l+2k (v)) we ob- 
tain decomposition (TTTb . □ 

Remark 4.2. The matrix coefficients of a polynomial F do only trivially depend 
on the choice of d. The matrix coefficients associated to d' > d will simply be the 
ones associated to d, enlarged by sufficiently many rows and columns of zeros. 
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Remark 4.3. From [2, Proposition 3.5], the polynomials P^it) are linear combi- 
nations of diagonal elements of the matrices with non negative coefficients. As 
a consequence, the matrix coefficients of any polynomial Pit) G K[i], are diagonal 
matrices. If Pit) = ^ fj.Pu(t), with all fk > 0, then the matrix coefficients Fk of 
P are also non negative, and, moreover, the top left corner of Fq equals fo- 

The following reformulation of Theorem 13. H is an analogue of the classical ex- 
pression of the linear programming bound (see e.g. (7J Chapter 9, Theorem 4]). 

Theorem 4.4. Let Eq be the matrix whose only non zero entry is the top left corner 
which contains 1. For a polynomial F(u, v, t) G let [Fo, . . . , Fj) be symmetric 
matrices such that 

d 

Fiu,v,t) = Y,{Fk,Y n k (u,v,t)). 

k=0 

Suppose the following conditions hold: 

(a) F k y for all < k < d. 

(b) F - f E y for some fo > 0. 

(c) F(u, v, t) < Ofor all (u, v, t) G A. 

(d) F(u, u, 1) < Bfor all u G [cos (f>, 1]. 

Then, for any code C in Cap(e, 4>) with minimal angular distance at least 9, 

card(C) < — . 

Jo 

Proof. The statement follows immediately from Theorem l3. 1 I because the matrices 
Gq = Fo/fo — Eq and G k = Fk/ fo for 1 < k < d are a feasible solution to the 
SDP GO} with M = B/fo - 1. 

We also give a direct proof, which has the additional feature to give information 
about the case when the obtained bound coincides with the size of a certain code. 
Let 

S := F(e ■ c,e ■ c r ,c- c r ). 

(c,c')ec 2 

We expand F in the Y^'s: 

S = E^' E ^(e-c,e-c', C . C ')). 

k=0 (c,c')eC 2 

On one hand, from the positivity property (© together with the fact that (A, B) > 
for any two positive semidefinite matrices A,B we. obtain 



S > (foE , E Yo{e-c,e-c,c-c)) 



(12) Mec* 



fo E ( y o) , ( e ' c ' e " c '' c " c ') = /ocard(C) 2 . 

(c,c')eC 2 
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On the other hand, if we split the sum S into diagonal terms belonging to pairs 
(c, c) and into cross terms belonging to pairs (c, d) with c ^ d, we obtain from 
condition (c) and (d) 

S = y^F(e ■ c,e ■ c, 1) + £\ c , )gc2 c ^ c , F(e ■ c, e • d, c ■ d) 
(13) ttc 

< Bcard(C) + 0, 

because (e ■ c, e ■ c, 1) G Ao and (e ■ c, e ■ d , c ■ d) G A if c 7^ c'. Now (fl2l and 
(TOT ) together give the inequality card(C) < B/fo. □ 

Remark 4.5. Life in LP method, the above proof gives additional information 
on the case of equality. Namely, if for a given code C and a given polynomial F, we 
have card(C) = B/ fo, the inequality (1131) must be an equality. So, F(u, v,t) = 
for all (u, v, t) running through the set of triples (e • c, e ■ d , c ■ d) with d and 
(c, d) G C 2 , and F(u, u, 1) = B for all u = e • c with c G C. 

Remark 4.6. In view of explicit computations, it is more convenient to remove the 
factor Xijfrom the coefficients ofYP, so that polynomials with rational coefficients 
have rational matrix coefficients. It changes the above defined F^ to congruence, 
hence does not affect the property to be positive semidefinite. These are the matrix 
coefficients we discuss about in the next two examples. 

Example 1. (d = 1) 

We consider the polynomial F = t — cos 6 — uv + cos 2 <p. The matrices of the 
decomposition (fTTT) are: Fq = (g J) with a = cos 2 (j) — cos 6 and F\ = (1). 
Condition (a) of Theorem 14.41 is fulfilled if a > 0. Condition (b) holds for /o = a. 
Obviously (c) holds if cos0 > and B = 1 — cos# because F(u, u, 1) = 1 — 
cos 9 — u 2 + cos 2 4>. We obtain: 

If cos 4>>0 and cos 6 < cos 2 4>, then A(n, 0,<p) < 



cos 2 6 — cos I 



It is worth to point out that the polynomial G = (t — cos 9) — cos 4>(u + v — 
2 cos 4>) leads to exactly the same bound. This time Fq = [ ~ l c ) with c = cos 4>, 
fo = a,B = l — cos 9. 

The above bound is already proved in [3, Theorem 5.2]. Indeed with the no- 
tations of (3l, let w (9, 4>) be defined by cos w(9,<j>) = (cos 9 — cos 2 (ft) / (sin 2 (ft); 
we have just proved that the Rankin bound for A(n — l,w(9, (ft)) is also a bound 
for A(n, 9, (ft). More generally, LP bounds for A(n — l,w(9, (ft)) are also bounds 
for A{n,9,(ft): Let f(x) be a polynomial of degree d that realizes an LP bound 
on S n ~ 2 for the angle w(9,4>). We can take polynomial approximations of the 
function 

F (^,^(( 1 -^)( 1 -^ ) )^ / ( ((i _j- p2))1/2 ) 

obtained by the truncated developments of the powers ( ( 1 — n 2 ) ( 1 — 1> 2 ) ) fc ^ 2 around 
u = cos 4>, v = cos <p. 
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Example 2. (d = 2) 

We consider the polynomial F = (t + l)(t — cos 9) + a((u — cos 4>)(u — 1) + (v — 
cos 4>)(v — 1)) . The parameter a > will be chosen later to optimize the bound. 
Condition (c) is obviously fulfilled and condition (d) holds with B = 2(1 — cos 9). 
The polynomial (t + l)(t — cos 9) has non negative coefficients if we expand it 
in terms of the basis PJ} (t) whenever cos 9 < 1/n. More precisely its constant 
coefficient equals — cos6>) while the two others are positive. So we only need 
to make sure that Fo is positive semidefinite. We find that: 

'2a(£ + cos0) + I -cos 9 -a(l + cos</>) o(l - £ 
F = | -a(l + cos^) (l-cos0) 



Let 



/ (o) := -a 2 ( (1 i +C ° S t )2 + (1 - -)) + 2a(i +co S ^) +-- cos 9. 



1 — cos 9 n J n n 

Then, an easy calculation shows that Fo ^ iff fo(a) > 0, and that Fo — /oFo >z 
iff /o < /o(o). The best bound is obtained when /o = /o(o) attains the maximal 
value 

,1 \ ( - + cos 6) 

(fo max = - " COS 9) + U ^ 2 ^ - 

v ra 7 / (l+cos t/i) 2 , _ j_ 

1— cosS n 



The final bound equals 



2(1 - cos i 



(fo) max 

and is valid as long as (/o)max > and (— + cos <fi\ > (this last condition holds 
because (/o)max must be attained at a positive a). 

It is worth noticing that the resulting bound is smaller than the LP bound for the 
entire sphere A(n, 9), obtained from the polynomial (t + l)(t — cos 9), which is 

2(1 - cos 9) 

For example, when cos <f> = cos 9 = 0, we recover the exact bound of 2n — 1 (see 
also 03). 

Remark 4.7. We can interpret the two examples treated above as follows: in both 
cases, we have perturbed the optimal polynomial for the LP method, respectively 
t — cos 9 and (t + 1) (t — cos 9), by a polynomial in the variables u, v, which affects 
the first matrix coefficient Fq and increases the value of the constant coefficient fo. 
However it seems difficult to generalize this approach. 

4.2. Orthogonality relations. In this subsection, we calculate the scalar product 
induced on R[u, v, t] by the natural scalar product on Po^S"™ -1 ) defined by ([Hi- 
Proposition 4.8. Let P € v, t] be a polynomial. We have 
1 f 



P(e ■ x, e ■ y, x ■ y)dw n (x)doj n (y) = I P(u,v,t)k(u,v,t)dudvdt 
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where 

K{U, V, I) = " " <..../,.. 



]fc( U , V ,t) = -n-^n-z r l _ u Z_ v 2_ { 2 + 2uvt )^ 



0, = {(u, v, t) : — 1 < u, v, t < 1, 

1 + 2twf -u 2 -v 2 -t 2 > 0}. 

Proof. If u = e ■ x and £ € 5 n_2 is denned by x = ue + (1 — u 2 ) 2(, we have 

du) n {x) = (1 — u 2 )^~ dudui n _i((). 
With y = ve + (1 — v 2 ) we have 

P(e ■ x,e-y,x ■ y)dio n (x) 

rl i _ 3 

/ P(u,v,t«; + ((l-u 2 )(l-?; 2 )) 2 C-0(l-^ 2 )~^^n-l(C) 

/•I /-l 

II / \ 9 \ n ~4 / n. re — 3 

= w n _ 2 y y P(«, v,t)(l - a ) 2 (1-u ) 2 dadn, 

where i := uv + ((1 — u 2 )(l — f 2 )) 2 <x With this change of variables having 

i 

Jacobian ((1 — u 2 )(l — w 2 )) 2 we obtain 
P(e -x,e-y,x ■ y)dw n {x) 

-i 

= uj n -2 P(u,v,t)(l - u 2 - v 2 - t 2 + 2uvt)^{\ - v 2 y^dudt, 
Jn(v) 

where 

Q{v) = {(it, t) : -1 < u,t < 1, 

1 + 2uvt -u 2 -v 2 -t 2 > 0}. 

Hence 

/ P(e ■ x,e ■ y,x ■ y)duj n (x)duj n (y) 

I \ / 9 9 9 \ 

= cj n _iw n _2 / P(u,v,t)(l — u —v —t -\-2uvt) 2 dudvdt 

□ 



Definition 4.9. Wfr/z f/ze notations of Proposition \4. 8\ the following expression de- 
fines a scalar product on M.[u, v, t]: 

(14) [F,G]= / F(u,v,t)G(u,v,t)k(u,v,t)dudvdt. 



From Proposition 14.81 it is the scalar product induced by the standard scalar /?ro<i- 
uct^onVo\(S n - 1 ). 
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The subspaces H?J are pairwise orthogonal. Consequently the matrix coeffi- 
cients of Yu{u, v, t) are pairwise orthogonal for [-, •]. Their norm is also easy to 
compute, and we obtain the following useful formulas: 

Proposition 4.10. 

(a) For all k, k' and all i, j, i' , f we have 



\i,j,k),(i',j',k') 



(b) For all symmetric matrices A, B and all k, k' we have 

(i6) i(A,n)AB,n>)] = 6k ' k :^[ B) - 

Proof. Obvious. □ 

4.3. Characterization of the positive definite polynomials. In view of Theorem 
14.41 we are concerned with the construction of polynomials satisfying condition 
(a). We prove in this subsection that this property is stable under multiplication. 
We start with a characterization of the set of polynomials satisfying (a) of Theorem 
031 



Definition 4.11. We say that the polynomial F(u, v, t) € M.[u, v, t] is positive def- 
inite if, for all finite C C S n ~ l , for all functions a : C — > M, 

(17) a(c)a(c)F(e-c,e-c',c-c')>0. 

(c, c ')ec 2 

The polynomials F(u, v, t) of the form 

d 

F(u,v,t)=J2(F k ,Y n k (u,v,t)) 

k=0 

with Fk y for all < k < d are positive definite in the above sense. Note that 
(fTTT ) is slightly stronger than the positivity property of the matrices Y£ proved in 
Theorem I2.2t the argument is essentially the same, as it follows from the equality 

£ a(c)a(c')Z%(c, c') = ( £ a(c)^(c)) ( ^ a(c)^(c)) * h 0. 
(c,c')ec 2 cec cgc 

We prove with next proposition that all positive definite polynomials in Rd arise in 
this way. 

Proposition 4.12. Let F(u, v, t) € Rd- Let (Fq, . . . , Fj) be symmetric matrices 
such that 

d 

F(u,v,t) = Y,(Fk,Y n k (u,v,t)}. 

k=0 

If F is positive definite, then F k y Ofor all < k < d. 
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Proof. Let F(x, y) = F(e ■ x, e ■ y, x ■ y). By compactness, F is positive definite 
if and only if for all / 6 Pol^™- 1 ), 

,f(x)f(y)F(x,y)du n (x)duj n (y) > 0. 

As a consequence, if Q(x) is any matrix, 

,{Q{x),Q{y))F{x,y)dw n {x)dw n {y) > 0. 



Let us fix k G {0, . . . , d} and let Abe a (<Z— k+1) x (cZ— symmetric, positive 
semidefinite matrix. Because of expression © of Z k , we can write (A, ZV;(x, y)*) 
in the form {Q(x), Q{y)). Hence, 



JA,Z^(x,y) t )F(x,y)dw n (x)(ko n (y) > 0. 



In terms of the scalar product [■ , ■] this is equivalent to 

[(A,Y n k ),F)>0. 

Since from (O [(A,Y%},F] = (h^y 1 (A, F k ), we have proved that (A,F k ) > 
for all A h 0, and so F fc ^ 0. □ 

Remark 4.13. This characterization of positive definite functions is in fact already 
proved in [4 , Section III] in a more general context: for compact spaces which are 
homogeneous under the action of their automorphism group, but not necessarily 
two-point homogeneous. The assumption that the group acts transitively is however 
not needed in the proof. 

Corollary 4.14. Let F,G £ Rd- If F and G are positive definite, then the product 
FG is also positive definite. 

Proof. From Proposition 14. 121 it suffices to consider the case F = (A, Y k ), G = 
{B,Y™), where A and B are positive semidefinite matrices. Again, we write 

(A,Z%(x,y¥) = (Q(x),Q(y)) and {B,Zf{x,yf) = (T(x),T(y)). With the 
formula 

(Q(x),Q(y))(T(x),T(y)) = (Q(x) ®T(x),Q(y) ®T(y)) 

we have 

a(c)a(c')F(c,c')G(c,c') 

(c,c')ec 2 

= £ a(c)a(c')(Q(c),Q( C '))(T(c),T(c')) 

(c,c')ec 2 

= Yl Wc)Q(c)®T(c),a(c')Q(c>)®T(c'))) 

(c, c ')ec 2 
= (Uc, U c ) > 
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with 

[/ c = 5>( C )Q(c)®T( C ). 

cec 

□ 

4.4. Reproducing kernels. We define the kernel Kd : R 3 x R 3 — > R by 

n 

(is) ^((u^.tj.c^y.o) :=x; /i r 1 (n(«.«.*).^(« / y.< / )>- 

fc=0 

Proposition 4.15. The kernel Kd is the reproducing kernel of the space Rd, i.e., 
for all F € Rd and all (u 1 , v',t') G R 3 we have 

(19) [K d (;(u',v',t')),F]=F(u',v',t'). 

Proof. It is straightforward from (fl6l) . 

□ 
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